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We investigate the two-pion interferometry in ultrarelativistic heavy ion collisions in the granular 
source model of quark-gluon plasma droplets. The pion transverse momentum spectra and HBT 
radii of the granular sources agree well with the experimental data of the v / ijviv = 200 GeV Au-Au 
and y/SNN = 2.76 TeV Pb-Pb most central collisions. In the granular source model the larger initial 
system breakup time may lead to the larger HBT radii -Rout, Rsidc, and i?i on g. However, the large 
droplet transverse expansion and limited average relative emitting time of particles in the granular 
source lead to small ratios of the transverse HBT radii -Rout/Rsidc- 



PACS numbers: 25.75.-q, 25.75.Gz 

Hanbury-Brown-Twiss (HBT) interferometry is a use- 
ful tool to probe the space-time geometry of the particle- 
emitting sources in high energy heavy ion collisions [lj- 
4] . The experimental results of the HBT measurements 
for the Au-Au collisions at the high energies of the Rel- 
ativistic Heavy Ion Collider (RHIC) indicate that it is 
hard to describe the source space-time by a simple evo- 
lution model [5 8]. HBT interferometry data provide 
strong constraints for the models of source space-time. 
Recently, the HBT measurement for the ^/sjvjv = 2.76 
TeV Pb-Pb most central collisions at the Large Hadron 
Collider (LHC) is performed Q. A consistent explana- 
tion to the HBT data of the LHC and RHIC experiments 
is required naturally for the source models, which will be 
helpful to understand the initial condition, source evolu- 
tion, and particle freeze-out in ultrarelativistic heavy ion 
collisions. 

In Refs. [IfljLLl]]; the granular source model of quark- 
gluon plasma (QGP) droplets [l2j | is developed to ex- 
plain the RHIC HBT data @, @T In tnis work we in- 
vestigate the two-pion HBT interferometry in ultrarela- 
tivistic heavy ion collisions in the granular source model 
of QGP droplets. Our results indicate that the granular 
source for the LHC Pb-Pb collisions may have the same 
initial droplet temperature and velocity formula as those 
for the RHIC Au-Au collisions, but a larger initial sys- 
tem breakup time. The consistent granular source model 
reproduces the pion transverse momentum spectra and 
HBT radii in the most central collisions of the RHIC 
BEE! and LHC [IIH experiments. 

In ultrarelativistic heavy ion collisions, the system at 
central rapidity may reach a local equilibrium at a very 
short time To, then fast expand in the beam direction 
(z-axis). Because of the initial fluctuation, the local- 



equilibrium system is not uniform in space |16l Il7| . It 
may form many tubes along the beam direction dur- 
ing the fast longitudinal expansion, and finally fragment 
into many droplets (see Fig.l of Ref. (lOII d ue to the 
"sausage" instability and surface tension [1(|. On the 
other hand, the rapidly increased bulk viscosity in the 
QGP near the phase transition may also leads to the sys- 
tem breakup |ia |. 

We assume that the system fragments and forms a 
granular source of many QGP droplets at a time to(> r o)- 
On the basis of the Bjorken hypothesis [19(, the longitu- 
dinal velocity and rapidity of the droplets at to are 

/, 1 to + Zo . . 

v dz = zo/to, r] = - log , (1) 

2 t - zq 

and the transverse velocity of the droplets may be ex- 
pressed as [111, [20| 
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where zq and po are the longitudinal and transverse co- 
ordinates of the droplet at the breakup time to, 7£j_o is 
the maximum transverse radius of the system at to- In 
Eq. ^ , the quantities ar and br are the magnitude and 
power parameters of the transverse velocity which will be 
determined by the data of particle transverse momentum 
spectra. 

In our model calculations, the initial droplet radius in 
droplet local frame satisfies a Gaussian distribution with 
standard deviation ad, and the initial droplet centers are 
assumed distributing within a cylinder along the beam 
direction by pjj G3 
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where 7?._l and MZ±_ are the initial transverse radius and 
shell parameter of the granular source [lfl,LLli . Because of 



the longitudinal boost-invariant in ultrarelativistic heavy 
ion collisions, we may obtain the initial coordinate zq of 
the droplet by the longitudinal distribution 
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The evolution of the granular source after to is the 
superposition of all the evolutions of the individual 
droplets, each of them is described by relativistic hydro- 
dynamics with an initial local energy density eo and the 
equation of state (EOS) of the entropy density with a 
cross over between the QGP and hadronic gas |2l| - |23| . 
The values of the EOS parameters used in the calcula- 
tions are taken as the same as in Ref. 11] . In order to 



include the pions emitted directly at hadronization and 
decayed from resonances later, we let the pions freeze-out 
within a wide temperature region with the probability 

E2 



the initial droplet temperature T = 200 MeV [ll| in the 
calculations. The corresponding initial energy density of 
the droplets, eo, is 2.24 GeV/fm 3 . The standard devi- 
ation ad for the droplet radius distribution is taken to 
be 2.5 fm. For the granular source in the RHIC colli- 
sions, To and the breakup time to are taken to be 0.8 and 
4.3 fm/c, and the source parameters 1Z± and A1Z± are 
taken to be 6.5 and 3.4 fm, respectively. For the gran- 
ular source in the LHC collisions, these parameters are 
taken to be To = 0.4 fm/c, to — 8.0 fm/c, 1Z± = 7.8 
fm, and A1Z± = 5.5 fm. By comparing the transverse 
momentum spectra of the granular sources with the data 
of y/IWN = 200 GeV Au-Au [H, [fl and y/s^ =2.76 
TeV Pb-Pb [15[ most central collisions, we fix the pa- 
rameters of the droplet transverse velocity ax = 0.65 and 
bx = 1-40 the same for both the two energy sources in Eq. 
@, where 1Z±q is taken to be 8.0 fm in the calculations. 



^ oc J dir e-^>/^ + (l- /dir ) 



x e 



-(T c -T)/AT dcc ^ (7; > t > 80 MeV) , (5) 



where T c is the transition temperature, /di r is a frac- 
tion parameter for the direct emission, ATd; r and ATd cc 
describe the widths of temperature for the direct and de- 
cayed pion emissions. They are taken to be T c = 170 
MeV, / dir = 0.85, T dir = 10 MeV, and T doc = 90 MeV 

m. 
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FIG. 1: (Color online) The pion transverse momentum dis- 
tribution of the granular source and the experimental data of 

and STAR 13) and 
I) most central colli- 



v^^ = 200 GeV Au-Au (PHENIX H3 
y/sTTiv = 2.76 TeV Pb-Pb (ALICE Hi 



In Fig. [TJ we show the pion transverse momentum 
spectra calculated with the granular source model and 
the exp erimental data of -^/s^vjv = 200 GeV Au-Au 
QJ, [pf and yfsWjj =2.76 TeV Pb-Pb QJ] most central 
collisions. Assuming that the systems fragment when the 
local energy density decreases at a certain value, we take 
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FIG. 2: The space-time distribution of pion-emitting points 
in the granular source model. 

Figure[2]shows the space-time distributions of the pion- 
emitting points in the granular source model. Because 
of the higher collision energy, the granular source for the 
LHC collisions forms (system energy density decreases to 
a certain value) at the larger to, and hence has a wider 
z distribution (see Eq. ((H)) and larger 1Zj_ value (wider 
po distribution, see Eq. ([3])). These lead to the wider 
distributions of the longitudinal and transverse source 
points for the LHC granular source than those for the 
RHIC granular source. 

In Fig. [3l we show the HBT radii i?out, -Rside and 
i?i ong in the "out" , "side" , and "long" directions [2J, [25| , 
and the chaotic parameter A of the granular sources as 
functions of the transverse pion pair momentum, kr = 
|Pit + P2t|/2- They are obtained by fitting the two- 
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FIG. 3: (Color online) (a) The two-pion interferometry results 
of the granular sources and the experimental HBT data for 
the most central collisions of ^snn — 200 GeV Au-Au 
and yi^T =2.76 TeV Pb-Pb 



pion correlation functions in different fcr regions with 
the formula 

C(g ut,g s idc,fflo„g) = 1 + Ae-'o-t^ut-^dc^dc-ftUflU, 

(6) 
in the longitudinal comoving system (LCMS) [4]. The 
HBT data of the most central collisions of -^/snn = 200 
GeV Au-Au (STAR [8]) and ^/sWn = 2.76 TeV Pb-Pb 
(ALICE Q) are also shown. In our calculations^we use 
the same cuts as in the experimental analyses [8|, Q , that 
is the pion rapidity is limited with \y\ < 0.5 for the gran- 
ular source for the RHIC collisions and the pion pseu- 
dorapidity satisfies \r)\ < 0.8 for the granular source for 
the LHC collisions, respectively. It can be seen that the 
HBT radii of the granular sources agree well with the 
experimental data. The wider distributions of the trans- 
verse and longitudinal source points for the LHC granu- 
lar source (see Fig. [5]) lead to the larger transverse and 
longitudinal HBT radii than those for the RHIC granu- 
lar source. In experimental HBT measurements, many 
effects, such as the Coulomb interaction between the fi- 
nal particles, particle missing-identification, source co- 
herence, etc., can influence the results of the A [irlj]- 
Because we do not consider these effects in our model, 
our A results are larger. 



The ratio of the transverse HBT radii -R ut/-Rsidc is 
related to the average relative emitting time of the two 
pions and the velocity of the source transverse expan- 
sion, At = (|£i — t 2 \) and v ST @, Q. Although the larger 
breakup time to for the LHC source leads to a larger 
particle-emitting time, the limited value of At for the 
granular sources [10Hl2l . [261 ] and the large droplet trans- 
verse velocity lead to the small R ou t/Rsidc results (see 
Fig. [3] (e) and (e')). In the upper panel of Fig. HI we 
show the average pion-emitting time t and the average 
relative pion-emitting time At for the RHIC and LHC 
granular sources. It can be seen that the average pion- 
emitting time for the LHC source is much larger than that 
for the RHIC source. However, there is only small differ- 
ence between the average relative pion-emitting times for 
the LHC and RHIC sources. In the lower panel of Fig. [H 
we show the average source transverse and longitudinal 
velocities as functions of the transverse pion pair momen- 
tum kr- The transverse velocity v ST increases with fcr, 
and the longitudinal velocity 7J SL decreases with fcy. For 
the collisions at the LHC energy the granular source has 
larger transverse and longitudinal velocities. 
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FIG. 4: (Color online) Upper panel: the average pion- 
emitting time and relative pion-emitting time for the gran- 
ular sources. Lower panel: the average source transverse and 
longitudinal velocities for the granular sources. 

In summary, we investigate the two-pion HBT inter- 
ferometry in ultrarelativistic heavy ion collisions in the 
granular source model of QGP droplets. The pion trans- 
verse momentum spectra and HBT radii of the granular 
sources agree well with the RHIC ^/snn = 200 GeV Au- 
Au and LHC ^/sWn = 2.76 TeV Pb-Pb most central col- 
lisions. Our results indicate that the granular source for 
the collisions at the LHC energy may have the same ini- 



tial droplet temperature and velocity formula as those for 
the collisions at the RHIC energy, but a larger initial sys- 
tem breakup time. The larger breakup time may lead to 
the larger longitudinal and transverse distributions of the 
source points, and hence lead to the larger longitudinal 
and transverse HBT radii. However, the average relative 
emitting time of the two pions is limited in the gran- 
ular source model. The limited relative pion-emitting 
time and larger droplet transverse expansion lead to the 
smaller ratio of the transverse HBT radii, i?out/-Rsidc- In 
our granular source model, there is correlations between 
the space-time coordinates and source velocities. The 
data of the particle transverse momentum spectra and 
HBT radii for the collisions at the RHIC and LHC en- 



ergies give strong constraints for the model parameters. 
The consistent explain of the granular source model to 
these experimental data will be helpful to understand the 
formation and evolution of the pion-emitting sources in 
ultrarclativistic heavy ion collisions. 
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